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Abstract 

Previous calculations of the thermal /^-function in a hot Yang-Mills gas at 
the one-loop level have exposed problems with the gauge dependence and 
with the sign, which is opposite to what one would expect for asymptotic 
freedom. We show that inclusion of higher-loop effects through a static 
Braaten-Pisarski resummation is necessary to consistently obtain the lead- 
ing term, but alters the results only quantitatively. The sign, in particular, 
remains the same. We also explore, by a crude parameterization, the effects 
a (non-perturbative) magnetic mass may have on these results. 
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1 Introduction 



The behaviour of the effective couphng constant as = 5'^/47r in QCD at high temper- 
ature or density has been discussed for a long time, starting with the renormahzation 
group equation (RGE) arguments of CoUins and Perry that decreases logarith- 
mically at high density due to asymptotic freedom. The idea of QCD as a gas of 
weakly interacting quarks and gluons at high T originates from this observation. It 
has later been questioned if it is correct to use the same decreasing as as the renor- 
malized coupling constant when computing general n-point functions with non-zero 
external momenta, as the simple scaling assumptions used in [0] do not hold when 
the external momenta introduce extra dimensionful parameters. The zero temperature 
RGE can only be expected to be useful when the typical momenta involved scale with 
the temperature 0, |^. Also, the argument in |I| assumes that there are no infra-red 
problems, which are now known to exist Therefore, several groups have explicitly 
calculated the T-dependence of the three-point function in QCD at high T and used a 
renormahzation group equation, with the temperature and the external momentum k 
as scale parameters 0], in order to derive the running of as with (T, k) [6 - 15]. Even 
if as was found to decrease logarithmically at high T it would not be enough to justify 
an ideal gas approximation of QCD since the typical expansion parameters asT / k and 
■s/alT / n still grow at high T. 

Various problems and ambiguities arose when calculating the thermal /5-function. 
It was recognized soon that the dependence of ^^(T, k) on T depends strongly on 
which vertex is chosen to renormalize a^, the other vertices being determined by Ward 
identities [P, |, |ll|, |12]. This prescription dependence exists also at T=0 when 



the momentum-space subtraction is used |]IB[. There was also some ambiguity in the 
results which depended on whether the imaginary time formalism (ITF) or a real time 
formalism was used [pT] , |12|, but this is now better understood |jl3|. Furthermore, for 
a given vertex the /3-function depends on the momentum prescription and differs, for 
example, when the collinear and the symmetric points are used, both at zero external 
energy. Another problem arose in that the result is also gauge fixing dependent 
which puts into serious question the usefulness of such an approach. It is, in fact, not 
at all surprising that the /^-function shows a gauge dependence when computed using 
the standard effective action ^ since it is not gauge invariant off-shell. Landsman 
therefore proposed [[1^] to use the Vilkovisky-DeWitt effective action [|17|, |I9| to 
calculate an explicitly gauge independent /9-function, though it would still depend on 
the external momentum prescription. Also a Wilson-loop approach has been used to 
compute a gauge invariant quark-ant iquark potential from which an effective coupling 
was defined [pOf- Such a definition is not directly related to the coupling considered 
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here. 

In this paper we follow the prescription of |J10|], and use the Vilkovisky-DeWitt ef- 
fective action to calculate the three-gluon vertex at the static and spatially symmetric 
point at momentum k and at temperature r for a SU{N) Yang-Mills gas. This ap- 
proach has recently been used in [14, 15| where the one-loop /5-function was calculated 
and the scaling in r and k, was analysed. The choice of the static renormalization point 
can be partially motivated by the fact that in the ITF it is only the zero Matsubara 
frequency modes that are soft and need resummation (see Section |^). It also elimi- 
nates the problem of choosing between analytic continuations (retarded/advanced or 
time/anti-time ordered) which have different soft contributions [pT| . 

Since the /^-function here is linearly related to the two-point function by a Ward 
identity one might naively expect that it would have a high temperature dependence 
like t^/k^. However, at the static point there is a cancellation and it is found that at 



one loop [IC, 14 



dg .,2l7r2r 

T-r- = 



dr 871^" 16 K ■ ^^"^ 
The leading linear contribution does not come from the hard part of the loop integral, 
responsible for a r^/zt^-term, but from soft loop momenta. Therefore, in the spirit of 
the Braaten-Pisarski resummation scheme it is not consistent to stop the calcu- 
lation at the one-loop order for soft internal momenta, but the resummed propagator 
and vertices must be used to get the complete leading contribution. The main purpose 
of this paper is to perform the resummed one-loop calculation and analyse the new 
result. We do not include any fermion contribution since it is subleading at high T. 



2 Perturbative expansion of the /^-function 

The RGE with the temperature and momentum (r, n) as parameters was first derived in 
[§] using the fact that the renormalized n-point functions are formally independent of 
the renormalization condition. We would like to relate this RGE to a direct calculation 
of the derivative of the three-gluon function. Let us first fix the notation and work in 
the Landau gauge in this section — it can be shown that results in this gauge, using 
the background field method, coincide with those results of the Vilkovisky-DeWitt 
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effective action [|T9l. The inverse of tlie full propagator is 



-l\ab 
) fiu 



abi 



n^,(p) = A^,n^(P) + 5^,n^(P) 



A 



p p 



B 



— R — ^ 



, V. = P'U. -U-PP.. 



(2) 



The four-velocity of the heat-bath is given by = (1,0,0,0). The Ward identities 
in the Vilkovisky-DeWitt effective action are particularly simple due to the off-shell 
gauge invariance, and the spatial part of the three-gluon vertex, for static and symmet- 
ric external momenta, can be related to the transverse part of the polarization tensor 
through 



rg'fc(r, n) = g [g.^ip - g), + g^.iq - r), + g^r - p)j] [ 1 + 



rabc 



(3) 



where P^ = {pq = 0, p) , = |p|^ etc., p^ = = r'^ = k^, and the dots stand for 
terms orthogonal to Pi, qj and r^. When po = we have 211-^ = —J^i^u using the 
Minkowski metric. The wavefunction and coupling constant renormalizations (AJJ — >■ 
zl^'^Afi^, g ^ ^^"^Qr) are performed at (r, k) so that 
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^■t{r,^) = gRiT,K)r'^\[g,^ip-q), + cyd.] + ... 



cabc 



■ (4) 



We now define an effective coupling constant g{T^ k) at another temperature T by 

rg',^(T, k)zI'\T, r) ^ g{T, K)r''^i^[g,,{p - q), + cycl.] + . . .| , (5) 

where zI^'^{T,t) = Z^^ (T') j Z\l'^ ir') is the rescaling of the field which is required in 
order to keep the normalization of the two-point function. In this way (^(r, k) measures 
the non-linearity of the theory. It is now straightforward to derive 



/3. = r 



dg{T,K) 
dr 



g ^dU^iT^K) 



dT 



(6) 



T=r 



using the renormalization condition in Eq.(HD. Similarly we find 



dgjr, k) 
dn, 



- I 
2 d\p 



d /n^(r,|p|) 



iPl 



(7) 



p|=K 
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For a perturbative calculation of j3r can be treated similarly) we need only the 
expression for n^(T, k) in the vicinity of the renormalization point, and we then use 
this in fixing the initial condition for the RGE. If the renormalization point is chosen 
appropriately, we could expect reliable results in some regime around (r, k) from a 
one-loop computation of Il^{T, k). 

We know that the leading hard thermal loops are non-local and if we want to include 
them through some resummation we also need non-local counter terms. Therefore, we 
write the action as 



C = -\tT{F') + \A,{ 



(8) 



and associate the first t^^^ with the "bare" propagator and consider the other one as 
a counter term. Then, we impose on the transverse and longitudinal part of t^^^ the 
one-loop hard thermal loop form 



7r^(po,p) 
7r^(Po,p) 



9 



1 - 



) 


"l-^ln 


Po+P 




2p 


Po-p _ 



9 [^,f^) 



6 



Po 

p2 



Po \ Po_ 
) 2p 



In 



Po+P 
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It is enough to introduce the momentum dependence in vr^'-^ from hard thermal loops 
in order to resum the leading g'^r^ contribution. The effective propagator, defined by 



* — l\ab 



(10) 



has an explicit r dependence which leads to a r-dependent UV divergence already at 
the one-loop level, since the iD^^^ contains contributions from an infinite sum of higher 
order diagrams. This r-dependence disappears when all diagrams to a given order are 



included and thus the problem can be pushed to arbitrarily high order by per- 
forming the renormalization to higher order. In our approach we only have to assume 
that this had been carried out at some r when renormalizing ]1^{t,k,). After taking 
the T derivative and the limit T ^ t everything is finite. We also note that vertices 
have hard thermal loop corrections so they should be treated in a similar manner by 
adding and subtracting the effective vertices in Eq.(^, but we do not need them in the 
approximation we are using (see Section H). 

Let us now analyse the perturbative calculation of using the renormalized La- 
grangian in Eq.(^. The high temperature expansion of Ii^{T,K) is an expansion 
in g'^T^/n'^ and g'^T/n. The g'^T'^/n'^ only comes from the hard thermal loops and 
for each such diagram there is a corresponding counter term g'^r'^ / with the opposite 
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sign generated by the last term in Eq. (H) . This is so because the counter term is chosen 
to be exactly the hard thermal loop contribution. The /3^-function is finally computed 
as the derivative of n^(T, n) with respect to T at T = r. If a diagram contains two or 
more hard thermal loops the leading g'^T'^/ and g'^r^/ terms factor out in such a 
way that after taking the derivative and T —>■ r they cancel. It then follows that in the 
perturbative expansion of dg{T, K)/dT at most one hard thermal loop contribute in each 
diagram and it is in fact an expansion in g'^r/ti only. The cancellations are identical 
to what was found for the 0^-model in |^ except that here we must use momentum 
dependent counter terms since the hard thermal loops are non-local. Also the usual 
way of simply using improved propagators to do loop calculations, without the RGE, 
does indeed resum the leading powers of g^T^ j . The difference is here that g itself 
is not a fixed zero temperature parameter but it is defined through the solution to the 
RGE. Therefore, the expansion is really in powers of g^ir^ K)TjK and its value depends 
on the solution of the temperature renormalization group equation. The possibility 
of performing a perturbative expansion at high r for fixed k depends on whether this 
combination increases or decreases at large r. 

Before doing any actual computation of the resummed /3-function it is interesting to 
discuss what kind of new terms one can expect and what their consequences would be. 
Let us therefore write 

= /5f + + = g'-{co + + c,^) , (11) 

assuming a high r expansion (r ^ k). The contribution from hard thermal loops is 
denoted by since it is generated by a longitudinal mass (see Section ^. In the 
expansion of in Eq.(|I]) there is no contribution from any hard thermal loop and 
we expect that the use of resummed propagators will supply the Ci term of relative 
order gr/n. The inclusion of /3j from a transverse "magnetic mass" of order g'^r, as 
discussed below, would generate the C2 term. We assume that the initial condition is 
given at a temperature tq ^ k while we still have (/^(tq, k)tq <^ k so that we can do a 
consistent perturbation expansion in g'^To/n. As r increases the solution to the RGE 
determines whether 5'^(r, k)t/k stays small enough for the perturbative expansion to 
remain valid. With the positive sign in Eq.(|^) for the bare one-loop /^-function the 
coupling constant diverges at some r implying that the expansion breaks down. If 
the sign had been negative the solution would go like gij^K) ~ (t/k)"^/^, implying 
that gr/n increases and has to be resummed while g^rjn goes to a constant and can 
be treated perturbatively if it is not too large. In the present case the bare one-loop 
calculation gives a divergent but resummation of grjK terms may change this. In 
particular, in Eq. (p!lD , if c\ is negative at large gT j k it dominates over the constant 
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term and the asymptotic form of g{T,K) is (r//t)^^/^. The factor gr/n still increases 
and needs to be resummed (as done with the momentum dependent counter terms in 
the temperature renormalization group equation) but the q'^t/k terms actually go to 
zero and the exact high r limit would be under control. We have found (see Section |^) 
that Ci is actually zero but there is a correction to the constant cq, though it does not 
change the sign of f3r for large qt/k. 

It is also interesting to see what happens if a magnetic mass is present; although 
such an effect is believed to be non-perturbative, we could crudely mimic such a term 
perturbatively by introducing some constant ttit ~ g'^r by hand as the position of the 
pole of the static transverse mode. Assuming that C2 is negative and dominates we find 
that g{T,K) ~ (t/k)-^^^. It may thus be inconsistent to assume that g'^r/n is large 
since it goes to a constant, and one would have to solve the renormalization group 
equation with the full (r, k) dependence. We found that C2 is gauge dependent and 
positive in the Landau gauge. Again, even if this correction would make jSr negative 
it is not consistent to separate out g'^r/n and subleading constants since they all go to 
constants. 

To summarize, a negative ci term would cure the problem of a divergent pertur- 
bative expansion of (3r, but the actual result shows only corrections to the cq and the 
sign remains positive leading to a divergent g{T, k) at some finite r. 



3 The resummed one-loop calculation 

To find the /3-function in the scheme described in Section |I| we need to compute the 
transverse part of the polarization tensor at one-loop using the effective propagators 
and vertices including hard thermal loop corrections. We shall perform the calculation 
in an arbitrary covariant background field gauge for comparison with other results and 
to see which terms are gauge independent, though the Vilkovisky-DeWitt approach 
prescribes the Landau gauge. The Feynman rules in the background gauge can be found 
in p6| and a one-loop calculation of the polarization tensor at finite T was performed 



m 



Z7[ . Let us start with /3^-function without resummation in a general covariant 
background gauge, parametrized with ^. The Landau gauge ^ = was considered in 



14, Iq] and the Feynman gauge ^ = 1 in |28|. We can extract the result for general ^ 



from the calculation in Furthermore, the leading t/k comes from the IR dominant 
part of the loop and is determined by the n = Matsubara frequency. For diagrams 
that are UV-convergent we can extract the linear t/k term by simply restricting the 
sum to ra = 0. Diagrams that are not UV-convergent have to be summed over all n. For 
the integrals we are dealing with it turns out that if the diagram is only logarithmically 
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divergent it is in fact enough to take the n = term to get the correct leading real 
part. There is an example in ||2^ where this does not work for the imaginary part. 
The expression needed for the one-loop polarization tensor in a general background 
field gauge, including the ghost contributions, is (we are using a different sign 



convention than |27]) 



Ii,,{K) = -g^NTY. I 



+ 



(27r)3 

P2Q2 



9,.^K{P) - 



iD^,{P) 



p2 



(12) 



where P + Q + K = Q and the bare three-point vertex is 
r,^^(P, Q, K) = g^p {P - Q)^ + g^(Q-K+ ^p\ 



g,a U-P-^Q 



(13) 



Calculating the transverse function Il^{po = 0,p^ = /t^) we find using the bare propa- 
gators in Eq. (|12D the following result for of Eq. ( pT]) : 



Stt^ Id k 



(14) 



For ^ = and ^ = 1 this coincides with [|T4| and respectively, and confirms the 
conjecture in |]15| that the difference between their result and that of is due to the 
gauge choice. 

The general one-loop calculation with effective propagators and vertices is difficult 
but in our case there are some simplifications. First we consider the external energy 
to be zero, and in the ITF only the n = internal modes need resummation since all 
other modes are hard. The effective propagators are thus only needed for zero energy 
and then they take the simple form 



46^ 



ab 



+ rrij' 



PiPj ■ 

p2 . 



1 XX, cPiP^ 



+ Till 



p. 



(15) 



where rriT and mi are transverse (magnetic) and longitudinal (electric) masses respec- 
tively. The longitudinal electric mass m| = ^g'^Nr'^ comes from the one-loop hard 
thermal loops, but the transverse magnetic mass is zero perturbatively. Here we 
try to estimate its effects by hand by inclusion of the term rriT ~ 0{g^T) in the prop- 
agator as a crude approximation to the true (non-perturbative) situation. Only the 
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pure gauge boson diagrams are effected by the resummation. The correction to e.g. 
the tadpole diagram is 



T 

2 4- 



„ *abcd j-^*cd^ „ abed r^cd i 



abed T)cd 



(16) 



where the star * denotes effective vertices and propagators. Each of the two terms 
in Eq.(|l^) is quadratically divergent and receives contribution from all Matsubara 
frequencies at high T. The difference however is only logarithmically divergent and to 
get the leading t/k term only the n = mode is needed. It then follows that 7*"^'^^ is 
only needed for zero external energy and then it reduces to the bare vertex. Similar 
simplifications can be done for the bubble diagram. The degree of divergence is reduced 
by two when subtracting the unresummed result and that is enough for using the n = 
approximation. We write the additional contribution to 11"^ from non-zero and 
like 

(17) 



5Ii'\zL. zt) = SlU'Xzl) + st'^U^izr) + (1 - O^'T^^'^' izT) , 
where z^ = mi/ n and Zt = mx/ 1^- The explicit expressions turn out to be 
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In this the following integrals have been used: 
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(19) 



The limit of for small z is given by 
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(20) 



while for large z we find 
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(21) 



It is worth noting that 5^11"'" is independent of the gauge parameter ^ and that it 
contains terms that are potentially dominant for large zl- However, it turns out that 
the leading terms cancel between the tadpole and the bubble diagram, and that 5^11"^ 
only contributes to Cq (and not to Ci) in Eq.(|ll]). When added to the bare one-loop 
result of Eq.(||) we find for ^ = 



-^N 

Stt^ k 16 



(22) 



The resummation has not changed the sign but the quantitative result to this order, 
showing that it was necessary to include these effects for a consistent calculation. The 
results of |T0| , p^ , |T5| are in this sense incomplete, but the general conclusions are 
correct since the sign remains unchanged. They could have been drastically changed 
if, for instance, the linear mi had come non-vanishing and negative (see Section 
When including rriT the large z limit gives 



Pr 



r /(e + 3)2 + 12 



16 



2 ^ 
7^ + — + 



K 8 16 



(23) 



where everything inside the [ ]-parentheses comes from the inclusion of rriT- In the 
Landau gauge Pr is positive even when including m-r. Even though it is possible 
to choose a gauge with a large enough ^ in order to stabilize the running of a.^ it 
seems rather artificial since we only can argue in favour of the ^ = gauge from the 
Vilkovisky-DeWitt approach. 

We have numerically solved the RGE in the high r limit at a fixed momentum scale 
K (i.e. neglecting the vacuum contribution and expanding in k/t) but using the exact 
dependence on and m^. The result is presented in Figs.(|l],|^). In the ^ = gauge 
(Fig. (la)) the coupling constant diverges at a finite temperature, just like without 
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Figure 1: The running coupling constant in a SU{2) Yang-Mills theory for the initial 
conditions g{l) = 0.5, 1.0, 1.5 and 2.0. Only the thermal contribution is included. The hard 
thermal loops are resummed and a magnetic mass uit = cg^T with c = 0.24, is included. 
The value of c is taken from the numerical simulations in p^/ . In the Landau gauge, ^ = 
(Fig. (la)), the coupling diverges at a Enite temperature and perturbation theory breaks down. 
In the i = 2 gauge (Fig. (lb)) the contribution from non-zero mT prevents the coupling from 
diverging. 
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Figure 2: The effects of including or excluding mi and ruT in the (3r-function. The 
parameters are the same as in Fig. with ^ = and g{l) = 1. 

resummation. If we choose ^ > 1, e.g. ^ = 2 as in Fig. (lb), the contribution from ttit 
changes the sign of /?r for large t/k. To see the effect of resummation in the ^ = 
gauge we have computed gij/ k) with and without the contribution from rriL and ttit 
(see Fig.(|^)). We find that the quahtative behaviour is not drastically changed by the 
resummation. The inclusion of has a tendency to increase the growth of g{T / 
while mT pushes the divergence to a higher temperature. 



4 Discussion 

The problems associated with a consistent calculation to this order of the /3-function 
concerning the "wrong sign" and the gauge dependence are reminiscent of the early 
one-loop bare calculations of the gluon damping constant at rest. Such calculations 
also gave the "wrong sign" , in that the modes were anti-damped, and the results were 
also gauge parameter dependent. The use of the Vilkovisky-DeWitt effective action to 
address the problem of gauge dependence in this case did not resolve the problem com- 
pletely, as the damping constant calculated in this formalism still had the wrong sign. 
Indeed, there were also arguments for choosing the background field Feynman gauge 
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^ = 1 as the "preferred" gauge based on the gauge invariant propagator of Cornwall 
p9|l , but this too gave the wrong sign for the bare one-loop damping constant and was 



quantitatively different than the Vilkovisky-DeWitt gauge ^ = 0. The resolution to 
these problems was later supplied by the Braaten-Pisarski resummation scheme 



where a gauge invariant and positive result is obtained to first order. Lessons from this 
could be drawn for this calculation of the /3-function. The results presented here in- 
dicate that higher-loop effects can change the result quantitatively, but the particular 
corrections considered here were not enough to resolve the problems of the wrong sign 
and of gauge dependence. This may mean that if the renormalization group equations 
in this form are to provide a useful tool a further resummation below the soft 0{gT) 
scale is needed to do a consistent calculation for the /3-function; as argued in |15|, the 
fact that the combination k/t appears means that a large temperature expansion is in 
a sense the same as probing the infrared behaviour. The need for such a further resum- 
mation in this context can also be seen when the simultaneous running of the coupling 
constant with temperature r and momentum scale k is investigated; from Eqs.(P,|^), we 
see that with the particular resummation investigated here the integrability condition 

^T-/^- = (24) 
dr an 

is not automatically satisfied. This particular problem could be cured in a somewhat 
ad hoc manner by having the T-derivative in Eq.(|^) act not only on the explicit T 
dependence arising from the Matsubara frequency sum but also on the implicit T de- 
pendence of the masses rriL and m-p. Doing so changes the results presented here only 
slightly quantitatively, however. One might thus expect that an improved resummation 
scheme, as well as addressing the problems of the sign of the /3-function and of gauge 
dependence, would also yield an integrable set of equations for (yf(r, k). The need for 
such a resummation beyond that of Braaten-Pisarski has also been recognized in the 
calculations of the damping rates of moving particles and of the production rates of 
soft photons ||3^. Whether such a scheme can be developed and can be used to give 



tractable results remains to be seen. 
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